Temperature and particle size dependence of equilibrium order parameter of FePt 
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First, second and third nearest neighbor pair mixing potentials for equiatomic FePt alloys were 
calculated from first principles by the Connolly-Williams method within the canonical cluster expan- 
sion formalism. It was demonstrated that the Connolly- Williams potentials (based on completely 
ordered states) and the Korringa-Kohn-Rostoker coherent potential approximation (KKR-CPA) po- 
tentials (based on completely disordered state) can be brought into very close correspondence to 
each other simply by increasing the magnitude of the strain-induced interactions added to the KKR- 
CPA potential. Using the mixing potentials obtained in this manner, the dependency of equilibrium 
Llo ordering on temperature was studied for bulk and for (approximately) spherical nanoparticles 
ranging in size from 2.5 to 6nm. The order parameter was calculated using Monte Carlo simulation 
and the analytical ring approximation. The calculated order- disorder temperature for bulk (1495- 
1514 K) was in relatively good agreement (4% error) with the experimental value (1572K). For 
nanoparticles of finite size, the (long range) order parameter changed continuously from unity to 
zero with increasing temperature. Rather than a discontinuity indicative of a phase transition, we 
obtained an inflection point in the order as a function of temperature. This inflection point occurred 
at a temperature below the bulk phase transition temperature and decreased as the particle size 
decreased. Our calculations predict that 3.5nm diameter particles in configurational equilibrium at 
600° C (a typical annealing temperature for promoting Llo ordering) have an Llo order parameter of 
approximately 0.84 (compared to a maximum possible value equal to unity). According to our inves- 
tigations, the experimental absence of (relatively) high Llo order in 3.5nm diameter nanoparticles 
annealed at 600° C or below is primarily a problem of kinetics rather than equilibrium. 



I. INTRODUCTION 

Self-assembled, monodispersed FePt nanoparticles are 
being intensively investigated for possible future applica- 
tion as an ultra-high density magnetic storage medium. 
In order to be useful as a storage medium, however, these 
particles, because of their extremely small volume, V, 
must have sufficiently high magnetic anisotropy, K Ul to 
withstand thermal fluctuations of the direction of mag- 
netization. This requires values of the thermal stability 
factor, (K u V)/(k-gT), of approximately 50. The parti- 
cles are usually produced by a "hot soap" process that 
yields a disordered f.c.c. solid solution alloy (e.g. Ref. 
1). Such particles are not useful for information storage 
in the as-made state because they are paramagnetic at 
room temperature due to their low magnetic anisotropy. 

Typically, the particles are annealed at a tempera- 
ture T ~ 600° C in order to induce an ordered Ll 
phased. The layered Llo phase* is known from stud- 
ies of bulk alloys to have an extremely high magnetic 
anisotropy (K u = 7 x 10 7 erg/cm 3 ). This value of mag- 
netic anisotropy would provide a sufficiently large ther- 
mal stability factor to make 3.5nm diameter particles vi- 
able for information storage. 

Unfortunately, it appears to be difficult to a achieve a 
high degree of long range atomic order in FePt nanopar- 
ticles with <4nm diameter by annealing at T <600°C 
(e.g. RefQ). One can consider two possible reasons for 
the fact that it has not been possible to obtain well or- 
dered small particles. First, the observed order may be 



low because the particle is not in its equilibrium state due 
to the slow kinetics at low temperatures. Alternatively, 
the equilibrium order itself may be low even at relatively 
low temperatures because of the small size of nanopar- 
ticles. The latter explanation was suggested in Ref. 0. 
There, the order-disorder phase transition temperature 
was estimated to decrease with decrease of particle size. 
For particle sizes less than 1.5 nm in diameter, the phase 
transition temperature was found to be below the typical 
annealing temperature (T ~ 600°C). Therefore, particles 
of diameter less than 1.5 nm were predicted to have no 
long range order in their equilibrium state at 600° C. This 
explanation is in qualitative agreement with experiment. 
The difference between the experimental (4nm) and the- 
oretical (1.5nm) critical size for disappearance of Llo or- 
der at 600° C was attributed to the neglect of nanoparticle 
surface effects. 

From our point of view, however, the results obtained 
in Ref. require verification because of the limitations of 
the theoretical models used in that study. Namely, the 
interatomic potentials in alloys usually are much more 
complicated and long-ranged than the nearest neigh- 
bor Lennard- Jones model that was used. In addition, 
the order-disorder phase transition temperature was es- 
timated in Ref. |2 by comparing the free energies of 
completely ordered and completely disordered states; 
whereas in reality, the ordered state approaches (with 
increasing temperature) the phase transition point be- 
ing not completely ordered. Also, the disordered state 
would be expected to approach the phase transition (with 
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decreasing temperature), not with a completely random 
atomic distribution but with an atomic distribution that 
has substantial short range order. Moreover, it is known- 5 
that there is no formal phase transition in a finite system. 

The aim of the present paper is to determine the rea- 
son for the low value observed experimentally for the Llo 
order parameter in small nanoparticles at typical anneal- 
ing temperatures through the use of theoretical models 
that do not have the above-described limitations. 6 

Presently, there are two main approaches for calculat- 
ing effective atomic interactions (called mixing potentials 
in this paper) that determine the configurational behav- 
ior of atoms in an alloy. One is the Connolly- Williams 
methodic, which is based on the first principles calcu- 
lation of the energies of a number of completely ordered 
structures. The second approach is based on the coherent 
potential approximation (CPA) for the completely disor- 
dered state (e.g. Ref. Il Of) . In the present paper we uti- 
lize and compare both approaches. The first-principles 
calculations that we performed in order to implement 
the Connolly- Williams method were carried out in the 
generalized gradient approximation to density-functional 
theory, using the VASP program package, 11 To obtain 
mixing potentials from the CPA we used data obtained 
within the Korringa-Kohn-Rostoker coherent potential 
approximation (KKR-CPA) 12 i 13 i 14 i 15 in Ref. El 

To study the temperature, concentration and size de- 
pendencies of equilibrium long-range order in FePt bulk 
and nanoparticles we used Monte Carlo simulation (uti- 
lizing the Metropolis algorithmic) and the analytical ring 
approximation 1 ^. 



Thus any configurational state of the two-component lat- 
tice gas model (i.e. any particular distribution of atoms 
over the sites) can be determined by C-^ variables only. 

The energy of any state can be expanded as (e.g. Sec. 
2 in Ref. W$) 

a,R Q i> Q 2Ri,R 2 

(3) 

where Eo, E^ and -E^'^ are the coefficients of the 

energy expansion. E^}'"^ can be considered as pair 
I \ Iv 

interactions between two atoms of types a\ and ai sit- 
uated at Ri and R2 sites, respectively. Expression © 
is usually called a cluster expansion^ In (J3J, the terms 
proportional to the powers of higher than second 
power are not taken into account. Terms corresponding 
to non-pair atomic interactions were found to be small in 
equiatomic FePt alloy (see discussion in Sec. IHIjl so that 
we excluded them from the begining. 

By the use of (J2J, one can exclude from J3J and 
get 

R i R!,R 2 

where V w (i = 0,1,2) we call mixing potentials (other 
names are "interchange energies" and "effective cluster 
interactions" ) : 

< ) =^R 1 -^R 1 +E^r,R 2 -^r,Rj. (5) 

R 2 



II. LATTICE GAS MODEL AND CLUSTER 
EXPANSION 

We consider an Fc-Pt alloy in the framework of the two- 
component A-B lattice gas model (e.g. A=Pt, B=Fe). In 
this modeliS, two types of atoms are distributed over the 
sites of a rigid crystal lattice. The atoms are allowed 
to be situated only at the crystal lattice sites and each 
site can be occupied by only one atom. The atoms inter- 
act through the lattice potentials and can exchange their 
positions according to Gibbs statistics. The lattice gas 
model is the most commonly used model for describing 
substitutional ordering in alloys (e.g. Refs. 

The configurational state of the lattice gas considered 
here can be defined by the set of configurational variables 

C R- 

£, a ( 1, if the site R is occupied by an a- type atom 
R 1 0, otherwise 

(1) 

where a =A,B and R is the site radius- vector. Because 
each site contains only one atom, we have 



t/(2) _ P A,A opA,B p B,B (f ,^ 

rti,rt 2 rti,rt 2 rti,rt 2 It, , . rt , 2 

is the energy of lattice gas with B atoms only. The 
knowledge of the mixing potentials is important because 
they determine the energetics of the lattice gas. A num- 
ber of methods have been elaborated for calculation of 
the mixing potentials. 

III. CONNOLLY- WILLIAMS METHOD 

Within the Connolly- Williams methodic, the mixing 
potentials are assumed to have the symmetry of the dis- 
ordered state, i.e. to be configurationally independent. 

Thus, Vj£ is independent of R and the pair mixing po- 

(2) 

tential Vj^ depends only on the difference Ri — R2: 

(2) 

The pair mixing potential Vj^ can be represented by its 

(2) 

values Vs for different coordination shells: 



C R + C K = 1 C K = 1 - C R- 



(2) 



(8) 
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where R s is the radius-vector connecting any given site 
with another site that belongs to the s-th coordination 
shell with respect to that given site. Generally, we define 
a coordination shell for a given site to be all of the sites 
that interact with a given site with the same value of the 
pair mixing potential. 

Within the Connolly- Williams method, the energies Ei 
(z = 1,2,..., N E ) of some N e structures (i.e. long-range 
ordered periodic atomic distributions sometimes called 
superstructures; for example Ll ) are calculated by first 
principles methods. Writing Eq. @ for each of those 
N £ structures, one can get the following system of linear 
equations 

Ei = v^ Cl + y. Si S vP ^ (9) 
* = 1,2,.. .J' 1 

where 

* = NfTW' (10) 

iVA and are the total numbers of A and B atoms in 
i-th structure, respectively, Sij are the structural coeffi- 
cients, and Cj is the concentration of A atoms: 

Ci = Nt/(N* + N»). (11) 

The energy V^ ' of a pure B crystal is also calculated 
by first principles methods. So the problem is to find 
N s + 1 unknown mixing potentials and V s {2) {s = 
1, 2, . . . , N s ) that give the best fit [through Eq. ©] to 
the energies £j (i — 1, 2, . . . , N e ) calculated by the first 
principles methods. 

Usually, the choice of the structures to be used is ran- 
dom with a predisposition towards the well-known ones 
that are more easily calculated. It is preferable that the 
structures be experimentally observable for the alloy un- 
der consideration and that they correspond to the ground 
state or at least a low energy state (e.g. Ref. l7||8|). How- 
ever, since all configurations are possible for Eqs. J3} and 
any convenient structure can be used for determina- 
tion of the mixing potentials. 

So, in contrast to the usual procedure, our criterion 
for the choice of the structures is only that they be lin- 
early independent. By linear independence we mean that 
the main determinant of the system @ i.e. the deter- 
minant of the matrix ||<Ji,j(cj — 5V,) + Sij\\ is not zero 
(i,j = l,2,...,N s ;Sis the Kronecker delta). Among the 
linearly independent structures we chose those with less 
than twenty atoms per unit cell in order to reduce the 
cost of the first principles calculations. Note that if all 
of the structures have the same composition (see below) , 
there is only one ground state (if it is not degenerate) and 
the usual procedure of using only (or mainly) the ground 
states within the Connolly- Williams method is inappli- 
cable. One may hypothesize that the use of non-ground 
states should expand the applicability of the resulting 
mixing potentials to higher temperatures. 



In the present paper we mainly (except Figs. I9I12|) 
study the Fe-Pt alloy close to equiatomic composition 
c = 0.5. Accordingly, we chose input structures to be 
of the same composition Cj = 0.5. The reason for do- 
ing so is the following. If we use structures of differ- 
ent compositions to find the mixing potentials, we imply 
that the mixing potentials are exactly the same for all 
compositions. However, the electronic structure of an al- 
loy system generally depends strongly on composition. 
(For example, at low temperatures the FePt alloy shows 
a strong dependence of magnetic properties on compo- 
sition: Fe3Pt, FePt and FePt3 are ferromagnetic, ferro- 
magnetic and antiferromagnetic, respectively.) The mix- 
ing potentials are mainly determined by the electronic 
structure of the alloy and, therefore, one may expect 
a compositional dependence of the mixing potentials as 
well. For example, such a dependence of pair mixing po- 
tential was found numerically within KKR-CPA in Ref. 
23. For such systems with strong compositional depen- 
dence of properties one may expect faster convergence 
of the Connolly- Williams method results with respect to 
increase of N e and N s when the considered structures are 
of the same composition. Such an approach corresponds 
to the canonical cluster expansion formalism 2 ^ and was 
applied e.g. in Ref. |2jJ 

Note that the Connolly- Williams method is usually ap- 
plied using structures of different compositions assuming 
compositional independence of mixing potentials (grand 
canonical cluster expansion formalism 24 ). In this ap- 
proach, the compositional dependence of alloy properties 
(e.g. the asymmetry of the phase diagram with respect 
to c = 0.5) manifests itself through the presence of non- 
pair mixing potentials. 26 In our approach the composi- 
tion dependence of alloy properties arises partly from the 
composition dependence of the mixing potentials. This 
approach is similar to that used in methods based on the 
study of the disordered state (e.g. CPA). Therefore, a 
direct comparison of Connolly- Williams and CPA poten- 
tials is possible (see below Sec. Ivj) . 

The first principles calculations were performed within 
the generalized gradient approximation to density- 
functional theory, using the VASP program package with 
mainly default settings,— All calculations were spin po- 
larized. The effect of lattice vibrations was omitted. All 
structures were totally relaxed including shape and vol- 
ume relaxation of the unit cell as well as the individual 
displacements of atoms within the unit cell. The densi- 
ties of fc-points using the Monkhorst-Pack mesh in the 
corresponding full Brillouin zones were chosen to be sim- 
ilar for all considered structures and approximately equal 
to ^/a)^ ( a * s ^ e f- c - c - lattice parameter), but the con- 
vergence of the results was checked to verify that it was 
sufficient in all cases. For the calculation of [see Eq. 
(fTOfl ] for each structure we used the quantity (energy 
of pure Fe) calculated separately using the same param- 
eters (both unit cell and VASP parameters) as for the 
corresponding structure. We believe that this approach 
diminished the error of the difference Ei — V 1 - ^ (and cor- 
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respondingly of £;) compared to the error of Ei and 
themselves, due to the systematic cancelation of errors. 

The Ll structure was included among the structures 
considered. In this case, after atom position relaxation, 
we obtained 3.848A and 3.77lA for the a and c lattice 
parameters of the corresponding tetragonal lattice, re- 
spectively (c/a = 0.980). For comparison the experimen- 
tal values are 3.847A and 3.715A (c/a = 0.966) M In 
addition, our calculated results showed the Ll ferromag- 
netic structure to be more stable (i.e. has lower energy) 
than the antiferromagnetic one in accordance with ex- 
periment. We believe that this good correspondence be- 
tween theoretical and experimental results confirms the 
adequacy of our VASP first principles calculations. 

The results of the application of our Connolly- Williams 
method to FePt are presented in Figs. 11141 and Table. 
U In Fig. [3 in analogy with Ref. |23, we present the 
dependencies of cross- validation^ (CV) and least-squares 
fitting (LSF) errors 



(CV) 2 = (23) 



23 
i=l 



VASP 



(LSF) 2 = (23)" 



23 

E 



.CW n 2 



_CW\2 



(12) 



(13) 



as well as of the phase transition temperature on N s 



In Eqs. D2Jl-©: 



e 



VASP 



(i = 1,2, ...,23) correspond to 



the values obtained within the VASP code for 23 linearly 
independent structures (see Appendix lAl ; ef w are ob- 
tained for those structures using the mixing potentials 
obtained by the Connolly- Williams method at N £ = 23; 
ef Wi are obtained for those structures using the mixing 
potentials obtained by the Connolly- Williams method at 
N e = 22 when i-th structure is excluded from fitting. 
Such a CV error characterizes the predictive power of 
these Connolly- Williams potentials within the set of 23 
struct ureA 

In Fig. the convergence of the results obtained by 
the Connolly- Williams method at N s = 3 is verified with 
respect to an increase in the number, N e , of structures 
taken into account within the method. Up to 23 linearly 
independent structures (see Appendix were consid- 
ered. Note that to achieve the rapid convergence shown 
in Fig. [3]it was important to begin with structures having 
the highest symmetries. The use of linearly independent 
structures guarantees that the good convergence (start- 
ing from N e ~ 8) is not caused by a simple similarity of 
structures. 

In Fig. 2] we compare the values of £j — V^Cj [see 
Eqs. (|9ll0|l ] for each of the twenty three linearly inde- 
pendent structures (see Appendix lA"| . The figure shows 
the values calculated both directly from first principles 
and by the use of the mixing potentials obtained by the 
Connolly- Williams method at N s = 3, N E = 23. The 
absolute differences between those two energies charac- 
terize the accuracy of Connolly- Williams fitting in the 
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FIG. 1: The values of pair mixing potential Vy for nine co- 
ordination shells (s = 1, 2, . . . , 9) corresponding to the KKR- 
CPA method (KKR-CPA), KKR-CPA plus strain-induced 
interactions (KKR-CPA+SI) at L = -0.0294; -0.0925 (see 
Sec. E) as well as to the Connolly- Williams method (CW) at 
N s = 3, N E — 23 (see Sec. 11111 . The distance is measured in 
f.c.c. lattice parameter units. 



case of each structure. The average accuracy of 5.72% 
(8.63 meV) per one structure was obtained. In order to 
check the predictive power of these Connolly- Williams 
potentials outside the set of 23 structures, we also cal- 
culated the energies of five additional structures. This 
set of five structures includes all distinguishable (but not 
identical to the previous twenty three ones) equiatomic 
structures based on two-cubic unit cells in the f.c.c. crys- 
tal lattice (see Appendix lA"|l . The average accuracy of 
the Connolly- Williams potential for the five additional 
structures and for all 28 structures were obtained to be 
7.60% (10.52 meV) and 6.06% (8.97 meV) per one struc- 
ture, respectively. It is important that if we calculate 
the Connolly- Williams potentials using all 28 structures 
(N £ = 28) we obtain almost the same average accuracies 
of the energy fitting for the first 23, last five and all 28 
structures: 5.92% (8.53 meV), 7.79% (10.36 meV) and 
6.25% (8.85 meV) per one structure, respectively. Thus, 
one may conclude that the comparatively lower accuracy 
of fitting for the last five structures is not caused by not 
including them into the fitting set. Such a high predic- 
tive power confirms our choice of using only linearly inde- 
pendent (see above) structures in the Connolly- Williams 
method, because the last five structures in Fig. 0] are 
linearly dependent (but not identical) to the first twenty 
three structures. One may suppose that the use of ad- 
ditional linearly dependent structures in the Connolly- 
Williams method would simply impose extra weight fac- 
tors for some of the initial structures. 

The good convergence as a function of N e shown in 
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experiment 




FIG. 2: The values of cross-validation (CV) and least-squares 
fitting (LSF) errors (see Eqs. I12H - H13|| ') as well as of the phase 
transition temperature To (within the ring approximation—) 
calculated as a function of the number of pair mixing potential 
values N s taking into account within the Connolly- Williams 
method (N e = 23). Dashed line corresponds to the experi- 
mental value—. 
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FIG. 3: The values of pair mixing potential V"i 2 ' for three 
coordination shells (s = 1,2,3) obtained by the Connolly- 
Williams method at fixed N s — 3 but for different numbers 
N E of structures (all of equiatomic composition c = 0.5) taken 
into account within the method. 



(2) 

TABLE I: The values of pair mixing potential 14 for sev- 
enteen coordination shells (s = 1,2, ...,17) corresponding 
to the KKR-CPA method (KKR-CPA), strain-induced in- 
teractions (SI) at L — -0.0294, KKR-CPA method plus 
strain-induced interactions at L = -0.0294 (KKR-CPA+SI) 
(see Sec. |VJ as well as calculated by the Connolly- Williams 
method (CW) at N s = 3, N e = 23 (see Sec. |TTTJ|. See also 
Fig. The cartesian coordinates of vector R are given in 
a/2 units, where a is the f.c.c. lattice parameter. Potential 
values are in meV units. 



Shell 


R 


|R] la 


KKR-CPA 


SI 


KKR-CPA+SI 


CW 


1 


110 


0.707 


185.855 


-9.879 


175.976 


87.69 


2 


200 


1.000 


-2.997 


-2.484 


-5.481 


-39.46 


3 


211 


1.225 


12.413 


0.326 


12.739 


15.85 


4 


220 


1.414 


-5.714 


1.173 


-4.541 




5 


310 


1.581 


1.727 


-0.724 


1.003 




6 


222 


1.732 


1.214 


0.091 


1.304 




7 


321 


1.871 


0.907 


0.133 


1.040 




8 


400 


2.000 


0.974 


-0.053 


0.921 




9 


411 


2.121 


0.147 


-0.188 


-0.041 






330 




-2.997 


0.438 


-2.560 




10 


420 


2.236 


-0.050 


-0.209 


-0.258 




11 


332 


2.345 


0.015 


0.030 


0.045 




12 


422 


2.449 


0.099 


0.060 


0.160 




13 


431 


2.550 


-0.075 


0.074 


-0.001 






510 




0.030 


-0.058 


-0.028 




14 


521 


2.739 


-0.006 


-0.073 


-0.079 




15 


440 


2.828 


-0.042 


0.163 


0.121 




16 


433 


2.915 


0.008 


0.019 


0.027 






530 




-0.008 


-0.058 


-0.066 




17 


442 


3.000 


-0.007 


0.023 


0.016 






600 




0.000 


-0.037 


-0.037 





Fig. |3 the high predictive power within and outside the 
set of 23 structures, the small CV and LSF errors of the 
Connolly- Williams fitting shown in Figs. 121 II and the 
position of the local minimum^ of CV error as a func- 
tion of the number of fitting parameters in Fig. Elsuggest 
that the neglect of nonpair mixing potentials and of pair 
mixing potentials outside three coordination shells (i.e. 
N s = 3) 30 as well as consideration of N £ = 23 linearly 
independent structures is sufficient to obtain adequate re- 
sults within the Connolly- Williams method for the case 
considered here. The small difference between the experi- 
mental bulk order-disorder phase transition temperature 
and that obtained using these mixing potentials^ (see 
Fig. [21 and below Tab. If HI in Sec. IVII|) provides indirect 
evidence of adequacy of the potentials for equiatomic Fe- 
Pt. In Fig. [T]and Tab. fflthe pair mixing potential values 
obtained by the Connolly- Williams method at N s = 3, 
N E = 23 (four fitting parameters) are shown. Uniform 
weight factors were used in the fitting. We obtained also 
V"W = 1.394 eV. 

Independent calculations within the KKR-CPA and 
the lattice-statics method also suggest that the atomic 
interaction tail outside the fourth coordination shell is 
weak (see below Sec. EJ, despite the "k — > non- 
analyticity"^* 3 ^* 3 ^ of the Fourier transform of the corre- 
sponding pair mixing potential. 
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first principles 

Connolly- Williams (7V=3, iV=23) 




5 10 15 20 25 30 
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FIG. 4: The values of e, ~ V (1) d [see Eqs. j£E3] for twenty 
eight equiatomic (c, = 0.5) structures calculated both directly 
from first principles and by the use of the mixing poten- 
tials obtained by the Connolly- Williams method (at N s = 3, 
N £ — 23). The differences between those two energies char- 
acterize the accuracy of Connolly- Williams fitting in the case 
of each structure. The first twenty three structures (see Ap- 
pendix^l are linearly independent and they (and only they) 
were used in the Connolly- Williams method. The last five 
structures (right side of the vertical dashed line) correspond 
to all distinguishable (but not identical to the previous twenty 
three ones) equiatomic structures based on two-cubic unit 
cells in the f.c.c. crystal lattice (see Appendix^. Those five 
structures were not used in the Connolly- Williams fitting, but 
they were considered in order to assess the predictive power 
of the obtained Connolly- Williams potentials. 



In Ref. Isil the Connolly- Williams method with struc- 
tures of different concentrations was applied to FePt. 
Consideration was not limited to the ground state struc- 
tures. One "hypothetical" phase LIi was considered. 
The interchange model was restricted to two coordination 
shells of atomic pair interactions neglecting non-pair in- 
teractions. Only six structures were considered making 
the number of unknown potential parameters equal to 
the number of input parameters (energies of structures) . 
The convergence of the results obtained by the Connolly- 
Williams method with respect to an increase in number 
of considered structures was not checked. Consideration 
of Fig. indicates that for such a small number of struc- 
tures convergence may not be achieved because of errors 
in the energy calculation. For all structures, the relax- 
ation of the unit cell volume was performed in a way 
that neglected possible unit cell shape relaxation (result- 
ing e.g. in its tetragonality) as well as neglecting the 
individual displacements of atoms within the unit cell. 
Using potentials calculated in this way, the phase tran- 
sition temperature for equiatomic FePt bulk alloy was 
found to be 2070K, which is much higher than the exper- 
imental one ~ 1573K. 

In Ref. HI the ratio V 2 {2) /v} 2) = -0.729 was derived 



through the use of first principles calculations. This is 

(2) (2) 

different from the value V 2 jV{ = —0.45 obtained by 
the Connolly- Williams method used in the present paper. 
It is difficult to determine the reason for this difference 
because details of the first principles calculations of Ref. 
l3fil are not reported. It is possible that the elastic con- 
tribution to the mixing potential neglected in Ref. 

is responsible for the difference. In Ref. |35j the value 
(2) 

of Vy = 96 meV (in our designation) was fit by repro- 
ducing the experimental bulk phase transition tempera- 

(2) 

ture. This implies a value of the quantity V 2 of = —70 
meV. For comparison, our Connolly- Williams potential 
values (see Tab. HJ are very similar for the first coor- 
dination shell but almost a factor of two smaller for the 
second one. In addition, for the third shell, the Connolly- 
Williams potential is nonzero. Of course, the experimen- 
tal bulk phase transition temperature is not reproduced 
perfectly by our set of Connolly- Williams mixing poten- 
tials which was obtained without fitting to experiment 
data (see below Tab. EH in Sec. fvTl^ . 

IV. Llo PARALLEL ANTIPHASE DOMAINS 

For the case of the Ll structure, there is a possibility 
of creating parallel antiphase domains, i.e. the regions 
where the same Llo structures are shifted with respect 
to one another by the distance a/2 along the z-axis per- 
pendicular to the Llo layers. The antiphase domains are 
usually created as a result of nucleation and growth of the 
Llo structure in different places as the initially disordered 
sample is cooled below the order-disorder phase transi- 
tion temperature. The antiphase boundaries (APB) be- 
tween the antiphase domains contribute to the total en- 
ergy making the existence of antiphase domains favorable 
or unfavorable depending on the sign of the antiphase 
boundary energy. 

It is easy to show that such parallel APB energies 
vanish in the case of only nearest neighbor interactions. 
Thus, first principles calculation of parallel APB energies 
is important because it allows one to test the adequacy 
of the nearest neighbor interaction model for FePt. In 
addition, knowledge of the parallel APB energy is help- 
ful for interpretation of Monte Carlo simulation results 
(see below Sec. IVI|) for the order parameter obtained by 
averaging in real space. 

In order to estimate the parallel APB energies, we cal- 
culated the values of E , E a , E 2a , i.e. the total energies 
of the pure Llo structure and of two structures composed 
of antiphase domains periodically repeated along the y 
direction with steps a and 2a, respectively (correspond- 
ing to the structures 1-3 in Appendix lA^l . The differences 
5E a and 5E 2a : 

5E a =E a - E , SE 2a = E 2a - E (14) 

can be considered as the total energies of a parallel APB 
in the corresponding two cases. The results of the cal- 
culations are presented in Tab. El Note that, besides 
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TABLE II: The energies 6E a and 5E 2a of parallel APBs in 
the Llo structure in the cases of a and 2a distance between 
parallel APBs, respectively. The values shown are per two 
atoms in one APB (in eV units) and were calculated by first 
principles using the VASP program packaged- In the relaxed 
case, we minimized the total energy with respect to the vari- 
ation of the size and shape of the unit cell as well as the local 
displacements of atoms within the unit cell. In the unrelaxed 
case, the atoms were assumed to occupy the sites of a rigid 
f.c.c. lattice with the lattice parameters assumed to be the 
same as in the relaxed case, neglecting all the local displace- 
ments of atoms within the unit cell. 

Unrelaxed Relaxed 



SE a 
8E 2a 
SE a - SE 2a 



0.223 
0.213 
9.80*10" 



0.213 
0.198 
1.59*10" 



the energies of the parallel APBs themselves, 8E a and 
6Ei a also contain some contribution from the interaction 
between parallel APBs. Such a contribution must de- 
crease with increasing distance between parallel APBs. 
The small difference between 6E a and SE 2a (see Tab. [Bj 
indicates that the interaction between parallel APBs in 
the two cases considered is small and short-ranged. It is 
somewhat larger in the relaxed case in accordance with 
the longer range of elastic interactions. The positiveness 
of SE a and SE2 a indicates that parallel APBs are not 
energetically favorable. The nonzero values of 5E a and 
bE-2a mean that the nearest neighbor interaction model is 
not adequate for a description of FePt alloy in accordance 
with results obtained in Sec. IIIII 



V. KKR-CPA AND STRAIN-INDUCED 
INTERACTIONS 

Elsewher e - 14 ! 15 , detailed procedures have been 
described for calculating effective atomic interac- 
tions (pair mixing potentials) within the KKR-CPA 
methodologyi2ii£. KKR-CPA-based effective atomic 
interactions are determined from the response of the 
electronic structure of the high-temperature, disordered 
alloy to small amplitude concentration waves. The 
results obtained within KKR-CPA in Ref. in the case 
of Feo.sPto.s alloy are presented in Fig. ^and Tab. [I] 

From Fig. []] and Tab. [Q one may conclude that the 
Connolly- Williams and KKR-CPA potentials (both cor- 
responding to equiatomic concentration c = 0.5) are sim- 
ilar but not identical. In both, only the three first coor- 
dination shells are substantially different from zero. The 
signs of the potential values for the first three coordina- 
tion shells are the same within both methods. Bowever, 
the absolute value of the Connolly- Williams potential for 
the first coordination shell is considerably smaller than 
that of the KKR-CPA. For the second shell, the situation 
is opposite. One may attribute the difference between the 
Connolly- Williams and KKR-CPA potentials to the dif- 
ferent structural and magnetic states used for the calcu- 



lations: completely ordered ferromagnetic for the former 
and completely disordered paramagnetic for the latter. 

For verification of the correspondence of both poten- 
tials to the real FePt alloy we calculated the bulk order- 
disorder phase transition temperature (within the ring 
approximation and by the Monte Carlo simulation) using 
both of the calculated mixing potentials. Then we com- 
pared the obtained values with the corresponding exper- 
imental one. The results are presented below in Tab. IIIII 
One can see that Connolly- Williams and KKR-CPA po- 
tentials underestimate and overestimate the experimen- 
tal phase transition temperature, respectively, by almost 
the same difference. 

Within the KKR-CPA calculation, the nuclei were con- 
strained to occupy ideal f.c.c. crystal lattice positions 
so the effect of inhomogeneous static atomic displace- 
ments is neglected. (The homogeneous ones are taken 
into account by minimizing the alloy's total energy with 
respect to the lattice parameter of the underlying lat- 
tice.) The inhomogeneous static atomic displacements 
are expected to be small for FePt alloy because of the 
small tetragonality obtained experimentally 27 and the- 
oretically (see e.g. above Sec. Illlf) . For verification, 
we calculated the corresponding strain-induced contri- 
bution to the pair mixing potential independently within 
the Khachaturyan 3 ^ semiphenomenological theory based 
on the microscopic Matsubara-Kanzaki-Krivoglaz 3336 - 3 ' 
lattice-statics method. This approach allows us to take 
into account the anisotropy and discrete (atomic) struc- 
ture of a crystal lattice in contrast with the macroscopic 
continuum theory of elasticity. 

To calculate the strain-induced interactions, we used 
the expression derived in Ref. |38| for the Fourier trans- 
form of the dynamical matrix. Platinum and iron were 
chosen as solvent and solute, respectively. The pair mix- 
ing potential is independent of this choice since it is sym- 
metrical with respect to interchange of A and B, see Eq. 
|JHJ| . For the calculations, we used the following numerical 
values: the mass of a Pt atom M = 32.40 x 10~ 26 kg; the 
lattice parameter ao = 3.92A; phonon frequencies: lon- 
gitudinal Ul = 27r x 5.80 and transverse ujt = 2tt x 3.84 
(all in Trad s -1 }^; elastic constants of Pt: C u = 3.580, 
da = 2.536, C 44 = 0.774 (all in 10 12 dyn/cm 2 ) 40 . The 
coefficient L of linear dependence of alloy lattice param- 
eter a on the Fe concentration c, 



a = a (l + L c), 



L - — — 

ao dc 



(15) 



c=0 



was calculated as L = — 0.0294 — . The concentration in- 
terval used for the calculation of L was chosen so that the 
Vegard rule is well fulfilled since it is a necessary condi- 
tion of applicability of the Matsubara-Kanzaki-Krivoglaz 
lattice-statics method 32 i 33 i 36 i 37 . 

The corresponding results for the strain-induced inter- 
actions are presented in Fig. ^and Tab. [H From them it 
follows that the strain-induced contribution to the pair 
mixing potentials is comparatively small as we expected. 
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FIG. 5: The dependencies of the Fourier transform Vj^ = 

J2 V$ exp(ikR) of pair mixing potential on the wave- 
R 

vector along the high-symmetry directions— within the cor- 
responding first Brillouin zone. The three cases correspond 
to (1) the pair mixing potentials obtained by the Connolly- 
Williams method (CW) at N s = 3, N e = 23 (see Sec. |TTTJ (2) 
pair mixing potentials obtained within the KKR-CPA and (3) 
KKR-CPA plus strain-induced interactions (KKR-CPA+SI) 
at L = -0.0294; -0.0925 (see Sec. EJ. The big bold point 

(2) 

and triangle correspond to V£ _ Q m two cases of taking into ac- 
count the strain-induced interactions, when the Fourier trans- 
form demonstrates nonanaliticity— . 



The addition of this contribution to the KKR-CPA po- 
tential decreases the corresponding phase transition tem- 
perature so that it is very close to the experimental one 
(see Tab. HT1. 

From Fig. □ and Tab. [Q one can also see that the 
addition of the strain-induced contribution to the KKR- 
CPA pair mixing potential moves it toward the Connolly- 
Williams potential values. To emphasize this effect, we 
increased the absolute value of the coefficient L from 
0.0294 to 0.0925 (thus artificially increasing the strain- 
induced interactions) in order to bring the KKR-CPA 
plus strain-induced potential as close to the Connolly- 
Williams one as possible. The same tendency is observed 
in the Fourier transforms of the corresponding pair mix- 
ing potentials (see Fig. [SJ. One may conclude that 
the mixing potential calculated for the disordered state 
(KKR-CPA) can be brought into close agreement with 
the mixing potential based on the completely ordered 
states (Connolly- Williams method) simply by increasing 
the magnitude of the strain-induced interactions (which 
is proportional to L 2 ). Perhaps such an increase of the 
strain-induced interactions roughly describes the depen- 
dence of mixing potential on the order parameter (con- 
figurational excitations of the electronic structur e 25 ! 43 ). 

Note from Tab. If HI that the phase transition temper- 
ature corresponding to the KKR-CPA potential plus the 



artificially increased strain-induced interactions is much 
lower than the experimental one. This may be explained 
by the roughness of our approach. Namely, by bringing 
the KKR-CPA to the Connolly- Williams one in the 
vicinity of G point (where the difference between them is 

(2) 

largest), we may cause too large an error in V£ in the 
vicinity of X point, which is mainly responsible for the 
phase transition because of minimum of there (see 
Fig. \jpk. 

From Fig. ^ an( i Tab. [I] it follows that the strain- 
induced potential is of short range despite the "k — > 
nonanaliticity" 7-32 33 of its Fourier transform. One may 
consider this fact as an additional proof of correctness of 
our application of the Connolly- Williams method in Sec. 
If III using the short range pair mixing potential. 

Below, in our Monte Carlo simulation of finite small 
systems, we exclusively used mixing potentials obtained 
by the Connolly- Williams method. We did so because the 
most interesting (for the aim of the present paper) tem- 
perature region corresponds to the highly ordered state. 
For this reason, we expect that the Connolly- Williams 
potentials which are based on the completely ordered 
state should work better in that region than the KKR- 
CPA potentials which are based on the completely disor- 
dered state. 

In Ref. the values of pair mixing potentials for the 
FePt bulk alloy were obtained within the tight binding 
linear orbital method. For c = 0.5, the values for the sec- 
ond and fourth coordination shell mixing potentials are 
close to the corresponding KKR-CPA ones from Ref. ^| 
but they are approximately a factor of two smaller for the 
first and third shells. Order-disorder transition temper- 
atures of 662K, 816K, and 435K were obtained by us for 
c = 0.25, c = 0.50, and c = 0.75, respectively, using these 
potentials within the ring approximation^ These transi- 
tion temperatures are much lower than the experimental 
ones (1087K, 1572K, and 1622K, respectively)^. The 
corresponding transition temperatures calculated in Ref. 
l45l are not so different from the experimental ones be- 
cause the mean-field approximation, which usually over- 
estimates the transition temperature^, was used. The 
magnitude of the discrepancy with experiment was un- 
derestimated also because an old value for the transition 
temperature at c = 0.75 from Ref. ^3 was used for com- 
parison. 



VI. MONTE CARLO SIMULATION 

For the Monte Carlo simulations we utilized the stan- 
dard Metropolis algorithm— We applied free and pe- 
riodic boundary conditions for the cases of spherical 
nanoparticlea^ and parallelepipeds, respectively. 

To obtain a rapid approach to the equilibrium state, 
which is the subject of this paper, the starting configu- 
ration for each temperature was chosen to be the com- 
pletely ordered one. We also used simplified kinetics in 
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our Monte Carlo simulation. Namely, we allowed any two 
randomly chosen atoms (not only nearest neighbors) to 
exchange their positions without an additional diffusion 
barrier. 

The absence of a monotonic change of mean values of 
all quantities being calculated in the Monte Carlo simu- 
lation was chosen as a criterion for the achievement of an 
equilibrium state. To fulfil such a criterion, it was gener- 
ally necessary to perform 1000 - 20000 Monte Carlo steps 
(i.e. the interchange of two atoms chosen at random) per 
site. After the equilibrium state was achieved, each calcu- 
lated quantity was averaged for subsequent Monte Carlo 
steps. When the oscillation amplitudes of such averages 
(considered as a function of the number of Monte Carlo 
steps performed) became less than 5% (of the absolute 
value of the corresponding average) during the last 10% 
of steps (from the total number of steps carried out at an 
equilibrium state) , the values of the averages calculated 
at the last Monte Carlo step were accepted as equilibrium 
ones. 

Wc define the equilibrium Llo order parameter rj as the 
statistical average of the maximum value among three ab- 
solute values of "directional" order parameters rj Xl rj y ,rj z : 

V=(v)nc> V = ™X-{\Vx\,\Vv\,\v»\} (16) 

where rji (i — x, y, z) is defined as the difference be- 
tween the Fe atom concentrations at odd and even crystal 
planes perpendicular to z-th direction, (. . is the sta- 
tistical average over the Monte Carlo steps. We chose this 
definition of 77 because of the equivalence by symmetry of 
the x, y, and z directions of Llo order. In addition, one 
can obtain an equivalent structure (at c — 0.5) by chang- 
ing the sign of which results in the exchange of Fe and 
Pt atoms producing a configuration that is equivalent by 
symmetry to the original one. During Monte Carlo sim- 
ulation, we observed fluctuations that cause transforma- 
tions between these equivalent states (i.e. fluctuations in 
the sign and direction of 77, see Figs. 16170 . This is in addi- 
tion to the usual statistical fluctuations within one such 
state. The Llo order parameter rj, defined in Eq. Q16(l 
takes into account the fluctuation induced transforma- 
tions between the equivalent states. Note that because 
of the above-discussed symmetry equivalence, we obtain 
(»7i)MC = for any i = x, y, z at any temperature, when 
the statistical average is taken over a sufficiently large 
number of Monte Carlo steps. Note that in Ref. [3^ an 
order parameter similar to Eq. (|16f) was introduced for 
the same reasons. 

Usually during Monte Carlo simulation we observed 
one of the "directional" order parameters to be much 
larger than other two (see e.g. Figs. I6I7|I . This means 
that the presence of perpendicular anti-phase domains is 
not favorable. Examination of the atomic arrangements 
within the Monte Carlo simulation showed no evidence 
of parallel anti-phase boundaries such as those discussed 
in Sec. IIVI The absence of anti-phase domains confirms 
the validity of our definition of the equilibrium Llo order 
parameter in Eq. (|l(j|) . 
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FIG. 6: The dependence of the FePt Llo order parameters 
on Monte Carlo steps during averaging in the equilibrium 
state during simulation of a spherical nanoparticle with di- 
ameter d = 3.0nm (see Tab. IIVI at temperature T — HOOK. 
rj x ,rj y ,rjz are the "directional" order parameters (see text), r\ 
is defined in Eq. HIGH . 
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FIG. 7: The same as in Fig. El but for the "directional" 
order parameters as the spherical nanoparticle with diameter 
d = 3.0nm (see Tab. HVfl approaches its equilibrium state at 
temperature T = 800K. 



Note that even with our simplified kinetics (see above 
Sec. IVI|) . we observed a slowing down problem in ap- 
proaching the equilibrium ordered state at low temper- 
atures. For example, in Fig. [7|one can see that, for a 
long time, the nanoparticle can be in a metastable state 
with order parameter (~0.6) lower than that (^0.8) in 
the equilibrium state. 
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FIG. 8: The temperature dependence of the FePt equilib- 
rium Llo order parameter rj obtained within the analytical 
ring approximation— for bulk ("bulk") as well as by Monte 
Carlo simulation for parallelepiped ("ppd") samples contain- 
ing N = 20 3 ,40 3 , and 60 3 atoms and the spherical ("sph") 
nanoparticle with d = 5.94nm diameter (see Tab. If VI at (or 
near) equiatomic concentration c = 0.5. 



VII. BULK PHASE TRANSITION 

To obtain an additional verification of the calculated 
values of mixing potential, we calculated the order- 
disorder phase transition temperature in the bulk FePt 
equiatomic alloy using these values. To do so, we cal- 
culated the temperature dependence of the FePt equi- 
librium Llo order parameter within the analytical ring 
approximation^ for bulk as well as by Monte Carlo 
simulation for parallelepiped samples containing N = 
20 3 ,40 3 , and 60 3 atoms and for spherical nanoparticles 
with 6.0nm diameter at (or near) equiatomic composition 
c = 0.5. The results are presented in Fig. [S] 

From Fig. [8] one may conclude the following. The 
ring approximation (which exactly corresponds to bulk, 
i.e. to an infinite sample) clearly shows a phase transi- 
tion temperature at which the order parameter rj drops 
to zero. Strictly speaking, in all the cases considered 
here of finite size samples (sphere and parallelepipeds) 
there is no phase transition in accordance with a general 
theorem^. The order parameter rj continuously changes 
from unity to zero with increasing temperature and in- 
stead of a phase transition we obtain an inflection point 
in the rj(T) curve. In the case of the parallelepiped with 
60 3 = 216000 atoms, the inflection point is very similar 
to the phase transition. Comparing the three curves for 
parallelepipeds of different sizes one can imagine that the 
inflection point transforms into a phase transition in the 
limit of infinite size. The position of the inflection point 
is often used to approximate the bulk phase transition 
point in Monte Carlo simulations of finite size samples. 



According to our calculations, the inflection point cor- 
responding to the largest parallelepiped simulation can 
be considered to be a good estimate for the phase tran- 
sition within the Monte carlo method. We plotted the 
temperature of the inflection point, T/v, as a function of 
N~ x / 2 and observed the limiting value as N~ x / 2 — > 0. 
This extrapolation of T/v to N — > oo hardly differs from 
the temperature of the inflection point for N = 60 3 . 

From Fig. 0it follows that Monte Carlo simulation for 
parallelepipeds with periodic boundary conditions ("no 
surface" ) is much more appropriate to the bulk behavior 
than Monte Carlo simulation for the spherical nanopar- 
ticle with free boundary conditions. The number of par- 
ticles N = 20 3 = 8000 in the smallest parallelepiped 
considered is very close to the number N — 8007 of par- 
ticles in the spherical nanoparticle (see below Tab. HVj l. 
This large difference in Monte Carlo simulation between 
parallelepiped and sphere is caused by the difference in 
shape and boundary conditions. 

TABLE III: The values of bulk order-disorder phase transition 
temperature measured experimentally (EXPERIMENT) 40 

(2) 

and obtained by the use of mixing potentials Vs 
corresponding to the KKR-CPA method (KKR-CPA), 
KKR-CPA method plus strain-induced interactions at 
L = -0.0294; -0.0925 (KKR-CPA+SI), Connolly- Williams 
method (CW), see Tab. Q] The theoretical values were calcu- 
lated within the ring approximation (ring) and by the Monte 
Carlo simulation (MC). 





To (K) 


EXPERIMENT 


1572 


CW 


1495(ring), 1514(MC) 


KKR-CPA 


1610(ring) 


KKR-CPA+SI (L = -0.0294) 


1552(ring) 


KKR-CPA+SI (L = -0.0925) 


997(ring) 



The temperatures, 1495 K and 1514 K, estimated 
for the phase transition within the analytical ring 
approximation 1 ^ and Monte Carlo simulation, respec- 
tively, are in close correspondence to the experimental 46 
1572 K. This 4% error demonstrates the adequacy of the 
mixing potential values at c = 0.5 calculated from first 
principles (i.e. without a fitting to experimental data). 

In order to test the applicability of these values over 
the entire concentration interval, we calculated the order- 
disorder phase transition temperature as a function of 
concentration within the analytical ring approximation 18 
using the mixing potential calculated for c = 0.5. The 
corresponding results together with experimental ones 
are presented in Fig. [5] (see also Fig. HOfl . From this 
figure one may conclude that the equiatomic mixing po- 
tential is only valid in the vicinity of c = 0.5. Outside this 
vicinity it is necessary to use different mixing potentials. 
The asymmetry of experimental phase diagram with re- 
spect to c = 0.5 implies that the mixing potentials should 
be concentration dependent and/or nonpair>2& According 
to our canonical cluster expansion formalism 24 , within 
the Connolly- Williams method one should calculate the 
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FIG. 9: Bulk configurational T — c phase diagram ob- 
tained within the analytical ring approximation— using mix- 
ing potentials calculated by the Connolly- Williams method 
at c = 0.5. Lines correspond to phase transitions from 
disorder to Llo or LI2 structures with decreasing temper- 
ature. Stars correspond to the phase transitions observed 
experimentally4-(see also Fig. 1101 
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FIG. 10: Experimental bulk configurational T — c phase 
diagram^ 



mixing potentials at any given concentration using struc- 
tures with the same (or close) compositions. Within the 
KKR-CPA, the concentration dependence of mixing po- 
tentials is obtained naturally. 



VIII. NANOPARTICLE SIMULATION 

The mixing potentials calculated in Sec. IIIII by the 
Connolly- Williams method were used for Monte Carlo 
simulations of the temperature dependence of the equi- 
librium LIq order parameter of spherical FePt nanoparti- 



TABLE IV: The characteristics of spherical nanoparticles 
used in Monte Carlo simulation: d - diameter; d/a - ratio of 
diameter to lattice parameter; N, iVp e , and iVp t - the num- 
bers of sites, iron and platinum atoms in the sample, respec- 
tively; c = N-p e /Np t - concentration; r\ and D - the value and 
dispersion (due to the thermodynamic fluctuations) of order 
parameter at 600° C, see Fig. 1111 cmax - the value of con- 
centration for which the concentrational dependence of Llo 
equilibrium order parameter achieves maximum, see Fig. 1121 
d(nm) d/a N Np e ATp t c rj D cmax 
5.94 15.63 8007 4007 4000 0.500 0.863 0.007 0.52 
3.50 9.21 1601 825 776 0.515 0.839 0.021 0.53 
3.00 7.89 1055 519 536 0.492 0.746 0.033 
2.50 6.58 603 299 304 0.496 0.704 0.045 0.56 




TEMPERATURE (K) 



FIG. 11: The temperature dependence of the FePt equilib- 
rium Llo order parameter rj obtained within the analytical 
ring approximation— for bulk ("bulk") using KKR-CPA plus 
strain-induced interactions at L = —0.0294 (CPA+SI) and 
Connolly- Williams mixing potential at N s = 3, N e — 23 
(CW) as well as by Monte Carlo simulation for the paral- 
lelepiped ("ppd") sample containing N = 60 3 atoms and for 
spherical ("sph") nanoparticles with 2.5nm, 3.0nm, 3.5nm, 
6.0nm diameters d at equiatomic (or near equiatomic, see 
Tab. IIV0 concentration c = 0.5 using Connolly- Williams mix- 
ing potential. 



cles with diameters of 5.94nm, 3.5nm, 3.0nm, and 2.5nm. 
In all cases the compositions are close to c = 0.5. The 
characteristics of those nanoparticles are listed in Tab. 

Note that the starting configuration for all cases con- 
sidered was chosen to be a completely ordered one with 
an Llo order parameter of unity. The deviations of con- 
centration from equiatomic c = 0.5 are due to the finite 
sizes and spherical shape of the nanoparticles. 

The results of the simulations are presented in Fig. 
ITTI For spherical nanoparticles, the tendencies in the 
behavior of order parameter are the same as in case of 
parallelepipeds (see above Sec. IVIIfl . Namely, in case of 
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finite size particles, instead of a phase transition point 
(see bulk curve) there is an inflection point. The smaller 
the particle size the smoother the rj(T) curve, the less or- 
der parameter at low temperatures (below and somewhat 
above the inflection point) and the higher the order pa- 
rameter at high temperatures (well above the inflection 
point). The inflection point temperature decreases with 
decrease of particle size. Note that in Ref. |35| the same 
tendencies were observed but the position of the heat ca- 
pacity maximum was interpreted as a phase transition 
point. Such an interpretation should be very approxi- 
mate in case of small finite systems because there is no 
formal phase transition in a finite system^. The same 
tendencies as in Fig. 1111 were observed in Monte Carlo 
simulation of CU3AU nanoparticles in Ref. 

Our calculations predict, for example, that 3.5nm di- 
ameter particles in configurational equilibrium at 600° C 
would have an order parameter rj — 0.84 compared to a 
maximum possible value of unity (see Tab. IIVI and Fig. 
lilt . Therefore, annealing at 600° C will not yield perfect 
order for 3.5nm diameter particles. Approximately 17% 
of the atoms will be on the wrong sublattices, even in 
equilibrium. The corresponding dispersion of rj due to 
the thermodynamic fluctuations is comparatively small 
(e.g. 2.5% for d = 3.5, T = 600°C, see Tab. HV|) . 

The " tail" of the order parameter at high temperatures 
in Fig. El is a consequence of thermodynamic fluctua- 
tions of the order parameter. The asymptotic behavior 
of the dispersion of such fluctuations can be estimated 
as ~ 1/VN (N is the number of atoms So the value 
of the high temperature "tail" of the order parameter in- 
creases with decreasing particle size. Note that we define 
the order parameter as in Eq. (|16fl . taking into account 
the fact that two configurations with different signs of 
the order parameter are physically identical. 

Experimentally, nanoparticles are created with some 
dispersion in the concentration. Therefore it is impor- 
tant to know the concentration dependence of the order 
parameter. We studied this dependency for the cases 
of bulk and spherical nanoparticles of different diame- 
ters at fixed temperature T = 600° C. The correspond- 
ing results are presented in Fig. El At each concen- 
tration we present the equilibrium order parameter value 
so that near concentrations c = 0.5 and c = 0.25; 0.75 
the curves correspond to Llo and LI2 order parameters, 
respectively. 

From Fig. El it follows that for small FePt nanopar- 
ticles we observe an asymmetry in the concentration 
dependence of the order parameter with respect to 
equiatomic c = 0.5. This asymmetry is a consequence 
of the finite size of nanoparticles. It should be empha- 
sized that for the calculation of order parameters in Fig. 
1121 at each concentration we used the same pair mix- 
ing potentials as at c = 0.5 that results in a symmetric 
order parameter for case of bulk. The concentration de- 
pendence of mixing potentials and/or appearance of non- 
pair mixing potentials at c 7^ 0.5 can make the bulk curve 
asymmetrical^ and make an additional contribution to 
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FIG. 12: The concentration dependence of the FePt equi- 
librium Llo and LI2 order parameters in bulk ("bulk") 
and spherical ("sph") nanoparticles with 2.5nm, 3.5nm, and 
5.94nm diameters d at fixed temperature T = 600° C. The 
results for bulk and spherical nanoparticles were obtained 
within the analytical ring approximation— and by Monte 
Carlo simulation, respectively, using mixing potentials cal- 
culated at c = 0.5 within the Connolly-Williams method. 

the asymmetry for finite particles observed in Fig. 1121 
One should consider the results in Fig. El to be more 
approximate the further they are from c = 0.5. 

IX. SURFACE EFFECTS 

In our study, we have used mixing potentials obtained 
for infinite bulk alloys and free boundary conditions to 
simulate the equilibrium configuration of finite size parti- 
cles. The presence of the surface will change the atomic 
potentials in the near-surface region compared to bulk 
potentials. This change may result in surface segrega- 
tion, leading, for example, to the tendency of Fe or Pt 
atoms to be preferably situated at the surface. In general, 
surface segregation should decrease the total Ll order 
parameter in nanoparticles, because the Llo order will be 
reduced at the surface. The surface effect on the order 
will depend on the ratio of depth of surface segregation 
to nanoparticle size. 

Analytical estimation of such surface effects is not 
straightforward and will be done elsewhere— Here we 
only describe the source of the asymmetry observed in 
Fig. El Namely, it is observed that in the presence 
of a free surface the Connolly- Williams bulk mixing po- 
tentials result in weak segregation of Fe atoms to the 
surface. This causes a depletion of Fe in the interior of 
the nanoparticle. Thus the effective concentration of Fe 
atoms within the nanoparticle but outside the surface 
segregation region will be less than one half. Therefore 
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it is necessary to add Fe atoms to the interior of the 
nanoparticle to achieve equiatomic composition outside 
the segregating surface and accordingly to achieve the 
maximum Llo order in that region. That is why the 
maximum of the Llo order parameter lies on the Fe-rich 
side of equiatomic composition in Fig. 1121 A similar ef- 
fect (but for the pseudo phase transition temperature in 
finite systems^) was observed in Ref. |35j for the case in 
which a strong segregation potential leading to the segre- 
gation of Pt to the surface was introduced during Monte 
Carlo simulation. That additional potential also results 
in a decrease of total order parameter in accordance with 
our description above. 

For nanoparticles chemically synthesized by methods 
such as the " hot soap process' — , the problem of the ef- 
fect of the surface on the interatomic exchange potentials 
is even more complicated because these nanoparticles are 
likely to have unknown atoms and molecules attached to 
their surfaces. It should also be noted that truncation of 
the mixing potentials at a surface is an uncontrolled ap- 
proximation. The fact that our simulations yield segrega- 
tion of Fe to the surface should not be taken as evidence 
that this would happen for a real nanoparticle. 



X. CONCLUSIONS AND DISCUSSION 

In the present paper, the pair mixing potentials (which 
determine the alloy configurational behavior) were calcu- 
lated up to the third coordination shell for equiatomic 
FePt alloys from first principles using the Connolly- 
Williams method (Sec. HTT)l . It was shown that the 
application of this method does not necessarily lead to 
significant non-pair interactions if all of the structures 
employed in the Connolly- Williams method correspond 
to the same composition (canonical cluster expansion 
formalism^) . The mixing potential was also calculated 
by adding the strain-induced part of the interactions to 
mixing potentials obtained from the KKR-CPA approach 
of Ref. (Sec. EJ. The application of the Connolly- 
Williams method at fixed composition allowed the di- 
rect comparison of these potentials with those obtained 
from the KKR-CPA (plus strain-induced part). The mix- 
ing potentials are shown to be similar with both giving 
values for the order-disorder phase transition very close 
to the experimental one. It was demonstrated that the 
Connolly- Williams potentials (based on completely or- 
dered states) and the KKR-CPA potentials (correspond- 
ing to completely disordered state) can be brought into 
very close correspondence to each other simply by in- 
creasing the magnitude of the strain-induced interactions 
added to KKR-CPA potential (Sec. 

It should be emphasized that it was not clear a pri- 
ori that either the Connolly- Williams method or the 
CPA approach could be trusted for a system as com- 
plicated as FePt. Magnetic alloys are problematic for 
the Connolly- Williams method because the mixing po- 
tentials are strongly dependent on the magnetic state of 



the system. Thus the chemical order and magnetic struc- 
ture are expected to be intimately related. On the other 
hand the CPA approach is based on concentration fluc- 
tuations in a hypothetical high temperature disordered 
state. We consider the fact that both methods are able 
to approximate the order-disorder temperature well with- 
out adjustable parameters and that the mixing potentials 
derived from the two approaches are rather similar to be 
a remarkable result. 

It was shown by first-principles calculations that paral- 
lel Llo antiphase boundaries are not energetically favor- 
able (Sec. IIV|> . This fact demonstrates the inapplicability 
of the nearest neighbor interaction model in accordance 
with Connolly- Williams and KKR-CPA results. The ab- 
sence of parallel Llo antiphase boundaries was confirmed 
by Monte Carlo simulation using the Connolly- Williams 
potential (Sees. IV 111 Villi) . 

In Sees. IVIIIVIlil the mixing potentials obtained by 
the Connolly- Williams method were used to investigate 
the dependence of equilibrium Llo ordering on tempera- 
ture for bulk and for spherical nanoparticles with diame- 
ters of 5.94nm, 3.5nm, 3.0nm and 2.5nm. These calcula- 
tions used both Monte Carlo simulation and the analyt- 
ical ring approximation. The calculated order-disorder 
temperature for bulk (1495-1514 K) was in relatively 
good agreement (4% error) with the experimental value 
(1572K). For nanoparticles of finite size, the (long range) 
order parameter changed continuously from unity to fi- 
nite asymptotic values with increasing temperature. The 
nonzero asymptotic values are a consequence of thermo- 
dynamic fluctuations of the order parameter and are pro- 
portional to 1/ v^/V where N is the number of atoms in a 
sample. Rather than a discontinuity indicative of a phase 
transition we obtained an inflection point (a precursor of 
a phase transition at large size) in the Llo order as a 
function of temperature. This inflection point occurred 
at a temperature below the bulk phase transition temper- 
ature and which decreased as the particle size decreased. 

According to our investigations, the experimental ab- 
sence of (relatively) high order in nanoparticles below 
600°C is primarily a kinetic problem rather than an equi- 
librium one. For example, our calculations predict, that 
3.5nm diameter particles in configurational equilibrium 
at 600° C (a typical annealing temperature for promoting 
Llo ordering) have an Llo order parameter of 0.84 (com- 
pared to a maximum possible value equal to unity). It 
should be noted that to rapidly obtain the correct equi- 
librium state, we used simplified kinetics in our Monte 
Carlo simulation. Namely, we allowed any two randomly 
chosen atoms to exchange their positions without an ad- 
ditional diffusion barrier. In a real alloy, the main mecha- 
nism of atomic diffusion is much slower because it consists 
of the exchange of the positions of atoms and neighbor- 
ing vacancies through energy barriers. Moreover, at each 
temperature we started the simulation from the com- 
pletely ordered state, whereas the actual nanoparticles 
are initially prepared in disordered state and transfor- 
mation from the disordered to the ordered state may be 
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much slower than the reverse one, especially at low tem- 
peratures. Nevertheless, even with our simplified kinet- 
ics, we observed a slowing down problem in approach- 
ing the equilibrium ordered state at low temperatures. 
In real nanoparticles this problem must be much worse. 
Kinetic acceleration methods such as irradiation and/or 
addition of other types of atoms may be useful in ac- 
celerating the formation of long range order. Harrell et 
elm^ have observed ordering on annealing at tempera- 
tures somewhat below 400°C by adding Ag or Au to the 
particles in the synthesis step. However, it appears that 
this ordering is also accompanied by agglomeration and 
sintering of the particles. It should be emphasized that 
all of the results presented in this paper correspond to 
equilibrium states and are, therefore, independent of the 
particular kinetic pathways that lead to these states. 

In Sec. IVIIII it was demonstrated that for small finite 
samples even for a composition independent, pair mix- 
ing potential one obtains an asymmetric order parameter 
as a function of composition with respect to equiatomic 
composition. In particular, for FePt nanoparticles, if we 
neglect the effect of the surface on the mixing potentials, 
the maximum of the Llo order parameter as a function of 
composition is shifted from c = 0.5 toward the Fe rich re- 
gion. The connection of such an asymmetry with surface 
segregation is discussed in Sec. IIXI 
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APPENDIX A: DEFINITION OF STRUCTURES 

The present appendix (Tab. Ivland Subsecs. 1-28) de- 
fines the structures used in this paper. Unit cell basis 
vectors are given in Cartesian coordinates in units of the 
fee lattice parameter. The Cartesian axes are directed 
along the edges of the fee cube. The third unit cell ba- 
sis vector is perpendicular to the first and second ones 
for all structures. All of the structures have equiatomic 
composition c = 0.5. The atom positions in the tables 
are before relaxation. In addition to numbering, we also 
include the structure designations that we used during 
our calculations. For the first, second and twenty fourth 
structures we put also the designations known in the lit- 



erature ("Ll ", "CH(40)" and "Z2", respectively) M The 
first 23 structures are linearly independent and were used 
in the Connolly- Williams method in Sec. IIIII The rest 
5 structures are linearly dependent (but not identical) to 
the first 23 ones. They were considered in order to check 
the predictive power of the Connolly- Williams potentials. 

TABLE V: The characteristics of the structures used in the 
paper, i is the structure number, iV u , c , is the total number 
of atoms in the unit cell of structure, Eu.c. is the energy 
of structure per its unit cell (as obtained from VASP calcu- 
lations) , V^ c _ is the energy (per unit cell) of the structure 
with the same unit cell but with Fe atoms only (see Sec. [nj , 
Sis are the structural coefficients (i=l,2...,28, s=l,2...,8) in 
the Connolly- Williams system of equations @. Note that 
the values of e, [see Eqs. iHJt- HLH l can be determined from 



this table as e, = (-Eu.c. — Ki. c. )/-^u. c . i where the val- 
ues on the right side of this equation correspond to the i-th 
structure. 
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1. 4cub(l), Ll 



Unit cell basis vectors Cartesian coordinates: 
(1.0,0.0,0.0), (0.0,1.0,0.0), (0.0,0.0,4.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.5,0.0), (0.0,0.0,1.0), 

(0.5,0.5,1.0), (0.0,0.0,2.0), (0.5,0.5,2.0), (0.0,0.0,3.0), 
(0.5,0.5,3.0); 

B: (0.5,0.0,0.5), (0.0,0.5,0.5), (0.5,0.0,1.5), 



15 



(0.0,0.5,1.5), (0.5,0.0,2.5), (0.0,0.5,2.5), (0.5,0.0,3.5), 
(0.0,0.5,3.5). 



2. 4cub(2), APB(l), CH(40) 

Unit cell basis vectors Cartesian coordinates: 
(1.0,0.0,0.0), (0.0,1.0,0.0), (0.0,0.0,4.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (0.5,0.5,1.0), 

(0.0,0.5,1.5), (0.0,0.0,2.0), (0.5,0.0,2.5), (0.5,0.5,3.0), 
(0.0,0.5,3.5); 

B: (0.5,0.5,0.0), (0.0,0.5,0.5), (0.0,0.0,1.0), 

(0.5,0.0,1.5), (0.5,0.5,2.0), (0.0,0.5,2.5), (0.0,0.0,3.0), 
(0.5,0.0,3.5). 



3. 4cub(3), APB(2) 

Unit cell basis vectors Cartesian coordinates: 
(1.0,0.0,0.0), (0.0,1.0,0.0), (0.0,0.0,4.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (0.0,0.0,1.0), 

(0.5,0.0,1.5), (0.5,0.5,2.0), (0.0,0.5,2.5), (0.5,0.5,3.0), 
(0.0,0.5,3.5); 

B: (0.5,0.5,0.0), (0.0,0.5,0.5), (0.5,0.5,1.0), 

(0.0,0.5,1.5), (0.0,0.0,2.0), (0.5,0.0,2.5), (0.0,0.0,3.0), 
(0.5,0.0,3.5). 



4. 4cub(4) 

Unit cell basis vectors Cartesian coordinates: 
(1.0,0.0,0.0), (0.0,1.0,0.0), (0.0,0.0,4.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.5,0.0), (0.5,0.0,0.5), 

(0.0,0.5,0.5), (0.0,0.0,1.0), (0.5,0.5,1.0), (0.5,0.0,1.5), 
(0.0,0.5,1.5); 

B: (0.0,0.0,2.0), (0.5,0.5,2.0), (0.5,0.0,2.5), 

(0.0,0.5,2.5), (0.0,0.0,3.0), (0.5,0.5,3.0), (0.5,0.0,3.5), 
(0.0,0.5,3.5). 



5. 4cub(5) 

Unit cell basis vectors Cartesian coordinates: 
(1.0,0.0,0.0), (0.0,1.0,0.0), (0.0,0.0,4.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.5,0.0), (0.5,0.0,0.5), 

(0.0,0.5,0.5), (0.0,0.0,1.0), (0.5,0.5,1.0), (0.5,0.0,1.5), 
(0.0,0.0,2.0); 

B: (0.0,0.5,1.5), (0.5,0.5,2.0), (0.5,0.0,2.5), 

(0.0,0.5,2.5), (0.0,0.0,3.0), (0.5,0.5,3.0), (0.5,0.0,3.5), 
(0.0,0.5,3.5). 



6. 4cub(6) 

Unit cell basis vectors Cartesian coordinates: 
(1.0,0.0,0.0), (0.0,1.0,0.0), (0.0,0.0,4.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.5,0.0), (0.5,0.0,0.5), 

(0.0,0.5,0.5), (0.0,0.0,1.0), (0.5,0.5,1.0), (0.5,0.0,1.5), 
(0.0,0.5,3.5); 

B: (0.0,0.5,1.5), (0.0,0.0,2.0), (0.5,0.5,2.0), 

(0.5,0.0,2.5), (0.0,0.5,2.5), (0.0,0.0,3.0), (0.5,0.5,3.0), 
(0.5,0.0,3.5). 



7. 4ppdl-2pl(l) 

Unit cell basis vectors Cartesian coordinates: 
(1.0,-1.0,0.0), (1.0,1.0,0.0), (0.0,0.0,2.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,-0.5,0.0), (0.5,0.5,0.0), 

(1.0,0.0,0.0), (0.5,0.0,0.5), (1.0,-0.5,0.5), (1.0,0.5,0.5), 
(0.0,0.0,1.0); 

B: (1.5,0.0,0.5), (0.5,-0.5,1.0), (0.5,0.5,1.0), 

(1.0,0.0,1.0), (0.5,0.0,1.5), (1.0,-0.5,1.5), (1.0,0.5,1.5), 
(1.5,0.0,1.5). 



8. 4ppdl-2pl(2) 

Unit cell basis vectors Cartesian coordinates: 
(1.0,-1.0,0.0), (1.0,1.0,0.0), (0.0,0.0,2.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,-0.5,0.0), (0.5,0.5,0.0), 

(1.0,0.0,0.0), (0.5,0.0,0.5), (1.0,-0.5,0.5), (1.0,0.5,0.5), 
(0.5,0.0,1.5); 

B: (1.5,0.0,0.5), (0.0,0.0,1.0), (0.5,-0.5,1.0), 

(0.5,0.5,1.0), (1.0,0.0,1.0), (1.0,-0.5,1.5), (1.0,0.5,1.5), 
(1.5,0.0,1.5). 



9. 9ppdl-lpl(l) 

Unit cell basis vectors Cartesian coordinates: 
(1.5,-1.5,0.0), (1.5,1.5,0.0), (0.0,0.0,1.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,-0.5,0.0), (1.0,-1.0,0.0), 

(0.5,0.5,0.0), (1.0,0.0,0.0), (1.5,-0.5,0.0), (1.0,1.0,0.0), 
(1.5,0.5,0.0), (0.5,0.0,0.5); 

B: (2.0,0.0,0.0), (1.0,-0.5,0.5), (1.5,-1.0,0.5), 

(1.0,0.5,0.5), (1.5,0.0,0.5), (2.0,-0.5,0.5), (1.5,1.0,0.5), 
(2.0,0.5,0.5), (2.5,0.0,0.5). 



10. 9ppdl-lpl(2) 

Unit cell basis vectors Cartesian coordinates: 
(1.5,-1.5,0.0), (1.5,1.5,0.0), (0.0,0.0,1.0). 
Atomic Cartesian coordinates 
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A: (0.0,0.0,0.0), (0.5,-0.5,0.0), (1.0,-1.0,0.0), 

(0.5,0.5,0.0), (1.0,0.0,0.0), (1.5,-0.5,0.0), (1.0,1.0,0.0), 
(1.5,0.5,0.0), (1.0,-0.5,0.5); 

B: (2.0,0.0,0.0), (0.5,0.0,0.5), (1.5,-1.0,0.5), 

(1.0,0.5,0.5), (1.5,0.0,0.5), (2.0,-0.5,0.5), (1.5,1.0,0.5), 
(2.0,0.5,0.5), (2.5,0.0,0.5). 



11. 9ppdl-lpl(3) 

Unit cell basis vectors Cartesian coordinates: 
(1.5,-1.5,0.0), (1.5,1.5,0.0), (0.0,0.0,1.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,-0.5,0.0), (1.0,-1.0,0.0), 
(0.5,0.5,0.0), (1.0,0.0,0.0), (1.5,-0.5,0.0), (1.0,1.0,0.0), 
(1.5,0.5,0.0), (1.5,0.0,0.5); 

B: (2.0,0.0,0.0), (0.5,0.0,0.5), (1.0,-0.5,0.5), (1.5,- 
1.0,0.5), (1.0,0.5,0.5), (2.0,-0.5,0.5), (1.5,1.0,0.5), 
(2.0,0.5,0.5), (2.5,0.0,0.5). 



12. Ippd2-lpl 



Unit cell basis vectors Cartesian coordinates: 
(1.5,-0.5,0.0), (0.5,1.5,0.0), (0.0,0.0,1.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (0.5,0.5,0.0), 
(1.0,0.0,0.0), (1.0,0.5,0.5); 

B: (0.5,1.0,0.5), (1.0,1.0,0.0), (1.5,0.0,0.5), 
(1.5,0.5,0.0), (1.5,1.0,0.5). 



13. 2ppd3a-lpl(l) 

Unit cell basis vectors Cartesian coordinates: 
(3.0,-1.0,0.0), (0.5,0.5,0.0), (0.0,0.0,1.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (1.0,0.0,0.0), 
(1.5,0.0,0.5); 

B: (1.5,-0.5,0.0), (2.0,-0.5,0.5), (2.5,-0.5,0.0), (3.0,- 
0.5,0.5). 



14. 2ppd3a-lpl(2) 



Unit cell basis vectors Cartesian coordinates: 
(3.0,-1.0,0.0), (0.5,0.5,0.0), (0.0,0.0,1.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (1.0,0.0,0.0), (2.0,- 
0.5,0.5); 

B: (1.5,0.0,0.5), (1.5,-0.5,0.0), (2.5,-0.5,0.0), (3.0,- 
0.5,0.5). 



15. Ippd4-lpl 



Unit cell basis vectors Cartesian coordinates: 
(2.0,-1.0,0.0), (0.5,0.5,0.0), (0.0,0.0,1.0). 



Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (1.0,0.0,0.0); 

B: (1.0,-0.5,0.5), (1.5,-0.5,0.0), (2.0,-0.5,0.5). 



16. 2ppd4a-lpl(l) 

Unit cell basis vectors Cartesian coordinates: 
(4.0,-2.0,0.0), (0.5,0.5,0.0), (0.0,0.0,1.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (1.0,0.0,0.0), (1.0,- 
0.5,0.5), (1.5,-0.5,0.0), (2.0,-0.5,0.5); 

B: (2.0,-1.0,0.0), (2.5,-1.0,0.5), (3.0,-1.0,0.0), (3.0,- 
1.5,0.5), (3.5,-1.5,0.0), (4.0,-1.5,0.5). 



17. 2ppd4a-lpl(2) 

Unit cell basis vectors Cartesian coordinates: 
(4.0,-2.0,0.0), (0.5,0.5,0.0), (0.0,0.0,1.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (1.0,0.0,0.0), (1.0,- 
0.5,0.5), (1.5,-0.5,0.0), (2.5,-1.0,0.5); 

B: (2.0,-0.5,0.5), (2.0,-1.0,0.0), (3.0,-1.0,0.0), (3.0,- 
1.5,0.5), (3.5,-1.5,0.0), (4.0,-1.5,0.5). 



18. 2ppd4b-lpl 



Unit cell basis vectors Cartesian coordinates: 
(2.0,-1.0,0.0), (1.0,1.0,0.0), (0.0,0.0,1.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (1.0,0.0,0.0), (1.0,- 
0.5,0.5), (1.5,-0.5,0.0), (2.0,-0.5,0.5); 

B: (0.5,0.5,0.0), (1.0,0.5,0.5), (1.5,0.5,0.0), 
(1.5,0.0,0.5), (2.0,0.0,0.0), (2.5,0.0,0.5). 



19. Ippd4-2pl 



Unit cell basis vectors Cartesian coordinates: 
(2.0,-1.0,0.0), (0.5,0.5,0.0), (0.0,0.0,2.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (1.0,0.0,0.0), (1.0,- 
0.5,0.5), (1.5,-0.5,0.0), (0.0,0.0,1.0); 

B: (2.0,-0.5,0.5), (0.5,0.0,1.5), (1.0,0.0,1.0), (1.0,- 
0.5,1.5), (1.5,-0.5,1.0), (2.0,-0.5,1.5). 



20. 2ppd3b-lpl 



Unit cell basis vectors Cartesian coordinates: 
(1.5,-0.5,0.0), (1.0,1.0,0.0), (0.0,0.0,1.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (1.0,0.0,0.0), 
(1.5,0.0,0.5); 

B: (0.5,0.5,0.0), (1.0,0.5,0.5), (1.5,0.5,0.0), 
(2.0,0.5,0.5). 
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21. Ippd2-2pl(l) 

Unit cell basis vectors Cartesian coordinates: 
(1.5,-0.5,0.0), (0.5,1.5,0.0), (0.0,0.0,2.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (0.5,0.5,0.0), 

(0.5,1.0,0.5), (1.0,0.0,0.0), (1.0,0.5,0.5), (1.0,1.0,0.0), 
(1.5,0.0,0.5), (1.5,0.5,0.0), (0.0,0.0,1.0); 

B: (1.5,1.0,0.5), (0.5,0.0,1.5), (0.5,0.5,1.0), 

(0.5,1.0,1.5), (1.0,0.0,1.0), (1.0,0.5,1.5), (1.0,1.0,1.0), 
(1.5,0.0,1.5), (1.5,0.5,1.0), (1.5,1.0,1.5). 



22. Ippd2-2pl(2) 

Unit cell basis vectors Cartesian coordinates: 
(1.5,-0.5,0.0), (0.5,1.5,0.0), (0.0,0.0,2.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (0.5,0.5,0.0), 

(0.5,1.0,0.5), (1.0,0.0,0.0), (1.0,0.5,0.5), (1.0,1.0,0.0), 
(1.5,0.0,0.5), (1.5,0.5,0.0), (0.5,0.0,1.5); 

B: (1.5,1.0,0.5), (0.0,0.0,1.0), (0.5,0.5,1.0), 

(0.5,1.0,1.5), (1.0,0.0,1.0), (1.0,0.5,1.5), (1.0,1.0,1.0), 
(1.5,0.0,1.5), (1.5,0.5,1.0), (1.5,1.0,1.5). 



B: (0.0,0.0,1.0), (0.5,0.5,1.0), (0.5,0.0,1.5), 
(0.0,0.5,1.5). 



25. 2cub(2) 



Unit cell basis vectors Cartesian coordinates: 
(1.0,0.0,0.0), (0.0,1.0,0.0), (0.0,0.0,2.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.5,0.0), (0.5,0.0,0.5), 
(0.5,0.0,1.5); 

B: (0.0,0.5,0.5), (0.0,0.0,1.0), (0.5,0.5,1.0), 
(0.0,0.5,1.5). 



26. 2cub(3) 



Unit cell basis vectors Cartesian coordinates: 
(1.0,0.0,0.0), (0.0,1.0,0.0), (0.0,0.0,2.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.5,0.0), (0.5,0.0,0.5), 
(0.0,0.5,1.5); 

B: (0.0,0.5,0.5), (0.0,0.0,1.0), (0.5,0.5,1.0), 
(0.5,0.0,1.5). 



23. 2ppd3-3c-2pl 

Unit cell basis vectors Cartesian coordinates: 
(3.5,-0.5,0.0), (0.5,0.5,0.0), (0.0,0.0,2.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (1.0,0.0,0.0), 

(1.5,0.0,0.5), (2.0,0.0,0.0), (2.5,0.0,0.5), (3.0,0.0,0.0), 
(0.0,0.0,1.0); 

B: (3.5,0.0,0.5), (0.5,0.0,1.5), (1.0,0.0,1.0), 

(1.5,0.0,1.5), (2.0,0.0,1.0), (2.5,0.0,1.5), (3.0,0.0,1.0), 
(3.5,0.0,1.5). 



24. 2cub(l), Z2 

Unit cell basis vectors Cartesian coordinates: 
(1.0,0.0,0.0), (0.0,1.0,0.0), (0.0,0.0,2.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.5,0.0), (0.5,0.0,0.5), 
(0.0,0.5,0.5); 



27. 2cub(4) 

Unit cell basis vectors Cartesian coordinates: 
(1.0,0.0,0.0), (0.0,1.0,0.0), (0.0,0.0,2.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.5,0.0), (0.5,0.0,0.5), 
(0.0,0.0,1.0); 

B: (0.0,0.5,0.5), (0.5,0.5,1.0), (0.5,0.0,1.5), 
(0.0,0.5,1.5). 



28. 2cub(5) 

Unit cell basis vectors Cartesian coordinates: 
(1.0,0.0,0.0), (0.0,1.0,0.0), (0.0,0.0,2.0). 
Atomic Cartesian coordinates 

A: (0.0,0.0,0.0), (0.5,0.0,0.5), (0.0,0.0,1.0), 
(0.0,0.5,1.5); 

B: (0.5,0.5,0.0), (0.0,0.5,0.5), (0.5,0.5,1.0), 
(0.5,0.0,1.5). 
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